6, 2008 21:54 WSPC/INSTRUCTION FILE kaz2005b 



International Journal of Modern Physics E 
© World Scientific Publishing Company 



ROTATION OF TETRAHEDRAL NUCLEI IN THE CRANKING 

MODEL 



N. SCHUNCKi'2, P. OLBRATOWSKI^, J. DUDEK^, and J. DOBACZEWSKI^ 

(1) Departamento de Fisica Teorica, Universidad Autonoma de Madrid, 
Cantoblanco, 28 049 Madrid, Spain 
(2) Institute of Theoretical Physics Warsaw University, Hoza 69, PL-00681 Warsaw, Poland 
(3) Institut de Recherches Subatomiques , IN2Ps-CNRS/Universite Louis Pasteur 
F-67037 Strasbourg Cedex 2, France 

Received (received date) 
Revised (revised date) 

The three-dimensional cranking model is used to investigate the microscopic aspects 
of the rotation of nuclei with the tetrahedral symmetry. Two classes of rotation axes 
are studied corresponding to two different discrete symmetries of the rotating hamil- 
tonian. Self-consistent Hartree-Fock-Bogoliubov calculations show that the tetrahedral 
minimum remains remarkably stable until the first single-particle crossing. 

1. Introduction 

The prediction that atomic nuclei can possess stable low-lying configurations with 
the tetrahedral symmetry has been confirmed by several independent calculations 
using a variety of nuclear mean-fields^'^'"^'^. However, no experimental evidence has 
been reported so far, which can be put down to the lack of understanding of the 
excitation mechanisms in such exotic systems. We investigate here the collective 
rotation, since the deep minima and large barriers reported e.g. in Ref.^ suggest 
tetrahedral nuclei may be amenable to sustain rotational bands. This article is a 
sequel to previous works where the structure and moments of inertia of tetrahedral 
collective bands were analyzed^, and the consequences of the collective rotation 
in terms of symmetries and quantum numbers were discussed^. To complement 
these studies, we wish to focus here on the stability of the tetrahedral minimum 
as a function of the rotational frequency, which is one of the pre-requisites to an 
eventual observation of rotational bands. 

Today there exists several articles in the literature about the nuclear tetrahedral 
symmetry (associated to the point group of symmetry T^) and we refer the reader 
interested in a general introduction to this subject to e.g. Ref.^. For the purpose 
of the present study, it suffices to recall that a nucleus with a non-axial 032 oc- 
tupole deformation (in the standard expansion of the nuclear radius on the basis 
of spherical harmonics, with deformation parameters ax^) has all the symmetries 
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of a regular tetrahedron. Although other realizations of the group T^, involving 
higher-order multipole terms with A > 7, also appear to generate stable minima in 
the potential energy landscape^, they will not be considered in this article. 

The cranking model was the subject of several comprehensive review articles and 
we refer to Refs.®'^° for a general discussion of its main features. We just recall here 
how the model should be adapted in the 3-dimensional case: the rotation is described 
with the help of 3 Lagrange multipliers lu = {ujx, ujy,ujz) which arc interpreted as the 
classical rotational frequencies along the x-, y- and z-axis of the body-fixed frame 
respectively. Equivalently, one can choose the spherical representation in which the 
rotational frequency vector Co is parameterized by (uj, (^). 

It is known that a collective rotation of a nucleus with prolate or oblate defor- 
mation takes place about an axis perpendicular to the symmetry axis. However, 
in the case of a nuclear shape with the Tj' symmetry, the quadrupole deforma- 
tion is equal to zero, and consequently no simple criterion exists to determine the 
sorientation of the rotation axis. 

A possible procedure to find optimum rotation axis consists in computing the 
total energy as a function of the orientation of the rotation axis, characterized by 
the two angles {9, (p), for different rotational frequencies w, as was done in Refs.^'^. 
One thus obtains two-dimensional maps whose minima signal the energetically- 
favored axes of rotation. These calculations were performed using a macroscopic- 
microscopic technique, in which the total energy is the sum of a liquid-drop con- 
tribution parameterized as in Ref.^^ and shell-correction extracted from a Woods- 
Saxon potential with the form defined in Ref.^^. The results^'^ suggest that at low 
rotational frequencies (up to w ^ 0.3 MeV/H), no particular ax;is of rotation is 
favored, while at higher frequencies, several well-defined minima emerge. However, 
these calculations assumed a fixed tetrahedral deformation and did not include 
pairing correlations: therefore, in the present study we want to see whether the 
tetrahedral minimum survives the increase of angular momentum. 

2. Stability of the tetrahedral minimum 

In order to study the stability of the tetrahedral minimum and to treat on the same 
footing the response of the nuclear system to rotation and pairing correlations, we 
performed the Hartree-Fock-Bogoliubov (HFB) calculations in the current imple- 
mentation of the code hfodd (v2.17k)^'' for ^^"Zr, which shows^ a pronounced 
tetrahedral minimum at / = 0. In our calculations we employed the SLy4^'* Skyrme 
force in the particle-hole channel and the density-dependent delta interaction in the 
particle-particle channel, whose intensity was fitted so as to reproduce the trend of 
experimental pairing gaps in the neutron-rich Zr isotopes, see Rcf.** for details. 

The calculations were performed for two families of rotation axes: (i) the axis 
passing through the middle of the edge of the tetrahedron, (ii) the axis passing 
through the tip of the tetrahedron. They will be referred hereafter as to edge and 
tip axes, respectively. These two axes are given a particular attention as they are 
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symmetry axes of the rotating tetrahedron; a 4-fold axis with inversion (edge) and 
a 3-fold axis (tip), i.e. these symmetries commute with the cranking hamiltonian. 

For the cdgx^ and tip axes, we performed the constrained HFB calculations. In 
the case of the edge axis we requested that the discrete antilinear T-simplexes, 

and Sy, be conserved. (We refer to Ref. for a comprehensive review of the 
discrete symmetries implemented in the code HFODD). This was equivalent to con- 
serving the z-signature symmetry, . They were supplemented by full symmetry- 
unconstrained calculations, which was the only option available for the rotation 
about the tip axis. In the latter case, we also used the option to readjust the orien- 
tation of the rotation axis self-consistently in the course of the iterations. We did not 
notice (i) any change in the orientation of uj fixed at the beginning of the iterations, 
at least up to the first single-particle crossing, (ii) any significative difference in en- 
ergy between the various orientations of lj and between the symmetry-constrained 
and symmetry-unconstrained case. This confirms a spherical-like behaviour of the 
rotating tetrahedron as claimed in Ref.^. 
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Fig. 1. Evolution of the deformation parameters in the tetrahedral minimum of ^^"Zr as function 
of the rotational frequency u> for the edge- (circles) and tip-axes (squares and diamonds) . For both 
classes of axes, the quadrupole deformation remains equal to 0. For the tip-axis, the theoretical 
value of ass in the exact limit is marked as the dashed line. 



With the standard choice of the spherical harmonics for the calculation of the 
multipole moments, the edge-axis coincides with the z-axis of the body- fixed frame^. 
The deformations extracted from the multipole moments can be plotted di- 
rectly, cf. Fig. 1. The curve with plain circles shows the "tetrahedral" deformation 
0:32. It is worth noticing that it is nearly constant as function of tu, thus proving 
the resilience of the tetrahedral symmetry to the rotation. 

In the case of the tip-axis, we chose to first perform a rotation of the nucleus 
at rest by the Euler angles (a, /?, 7) in such a way that the new tip-axis coincides 
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with the z-axis of the reference frame. This corresponds to an angle (3 such that 
cos^/3 =1/3, and a = 7 = 0. The multipole moments Qx^ must be transformed 
according to the general relation: 



where the -D*^/ (a, /3, 7) are the Wigner matrices. In our case, the initial deformation 
was characterized by Q32 7^ and all other moments were null. After the rotation, 
only QsQ and Qss arc different from zero, and they are related by: Q30/Q33 = 
— •\/T0/2. We report in Fig. 1 the equivalent aao and a^s deformations as function 
of the rotational frequency, together with the value of 0:33 in the exact symmetry 
limit. The little deviation from the exact symmetry case seems to increase slightly 
with Lu although remaining very small. 

In conclusion, we have investigated the rotation of tetrahedral nuclei using 
a fully-microscopic 3-dimensional cranking model. The results of self-consistent 
Skyrme Hartree-Fock Bogoliubov calculations show that the tetrahedral minimum 
remains as function of the rotational frequency and that practically no quadrupole 
polarization remains. The latter observation implies strongly hindered stretched- £'2 
transitions for the eventual rotational bands. 
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